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Abstract 

This paper studies the thermal properties of three-terminal mesoscopic dielectric systems in the nonlinear 
response regime at low temperature. For a symmetric three-terminal system, when the temperature is finitely 
different between the left and right thermal reservoir, the temperature of the central thermal reservoir is 
always higher than the averaging temperature of the others. This nonlinear thermal phenomenon is also 
observed for asymmetric three-terminal systems. At the end, a model of thermal rectification is presented. 
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With the development of the modern electronics, heat conduction, as the counterpart of electric 
conduction, has attracted much attention in recent years. Some research works yj, |2j, [j, |4|] focus 
on the universal quantum of thermal conductance which is predicted by Rego and Kirczenow Q5[] 
and has already been verified by experiment [6] . The quantum of thermal conductance indicates 
that the heat conduction is determined by the ballistic transmission of the acoustic phonon at low 
temperature. On the other hand, thermal rectification to control the heat flux [7, S. [loi [ill 12! 



temperati 



is very interesting. Experimental work to demonstrate the thermal rectification is reported 
recently [15]. 

Based on the above efforts, this paper further studies thermal rectification in ballistic heat con- 
duction. For heat conduction in a two-terminal mesoscopic dielectric system at low temperature, 
because transmission coefficients of ballistic phonons are independent on the temperature, there 

n 

should be no thermal rectification H12H . But, for electric conduction within a three-terminal bal- 
listic junction, previous investigations have indicated the nonlinear ballistic transport of electrons. 
That is, if voltages Vl and Vr are applied on the left and right branches of a symmetric three- 
terminal ballistic junction in push-pull fashion, with Vl = -Vr, the voltage at the central branch is 



always Negative 



jun ction in pus 



20, 



21, 



220 . This nonlinear 



second-harmonic generation, and logic function 11231. l24l 12511 . Motivated by the nonlinear electri 



property can be used for rectification, 



cal properties in three-terminal ballistic junctions, in this work, we study the thermal properties of 
three-terminal mesoscopic dielectric systems in the nonlinear response regime and try to propose 
a model of thermal rectification. The model works at low temperature in order to keep the ballistic 
transmission of phonons, similar with the model taken in Ref. |l4fl . 

The geometry of the symmetric three-terminal mesoscopic dielectric system is sketched in 
Fig. Q] Regions -L, -R and -C are left, right and central terminals, respectively. Region - J is the 
midsection. Assuming that the terminals are perfect and phonons coming from thermal reservoirs 
are not scattered within the terminals, the energy flux Q, from terminal i (i = L,R, C) flowing into 



the midsection J can be expressed as [2' 



2I 



27] 



2, = ^^ f [n(T i ,to)-n(T j ,to)]htx>T jiM (to) < ^, 



(1) 



where n(T h to) = [exp(hto/k B Ti) - l]" 1 is the Bose-Einstein distribution function of the phonons in 
the zth reservoir, T t is the equilibrium temperature of thermal reservoir i, oj im is the cutoff frequency 
of mode m in terminal i, t 7 - ; >(o>) = d{a>-a>j n )Tj^ nm and r yij , !m is the transmission coefficient from 
mode m of terminal i at frequency to across all the interface into the mode n of terminal j. 



Firstly we let the energy flux flows through the two terminals L and R and keep zero flux in 
terminal C (i.e. Qc = 0) and then calculate the temperature T c . In the linear response regime, 



the flux in terminal C can be written as [|26l 



m Qc = Zj=L,RG jC (T c - Tj) = 0, where G jC 



is the two-terminal thermal conductance from central terminal C to terminal j (= L or R). Due 
to the symmetry of the three-terminal mesoscopic dielectric system and the independence on the 
temperature of transmission coefficients, the thermal conductances from central terminal to left 
terminal and right terminal are equal, i.e. G LC = G RC . So the temperature of central thermal 
reservoir in the linear response regime is simply the averaging temperature of the left and the right 
thermal reservoirs, i.e. T c = (T L + T R )/2. 

When the system is not in the linear response regime and with the finite temperature difference 
of 2|Ar| between the left thermal reservoir and the right thermal reservoir, what is the temperature 
T c of the central thermal reservoir? To figure out it, we let T L = T + AT and T R = T - AT. 
By using the Taylor expansion of the Bose-Einstein distribution function n(T h co), when | A7"| is a 
small value, the temperature T c can be written as 

r c = r + ^(Ar) 2 + O[(Ar) 4 ], (2) 



where 



Here, trc,«i = Thc,m by the symmetry of the system and the independence on the temperature of 
the transmission coefficients. Thus, T c depends quadratically on AT for small \AT\. By using 
e x = 2Xo T\ x ^ f° r a U <^ > 0, we can easily obtain that 

l dn{T, oj) \ nu)exp(Hu)/k B T Q ) 
\ dT j To ~ k B T*[exp(tuo/k B T ) - l] 2 ' 
d 2 n(T, a>) \ fioj exp(ha)/k B T ) 



dT 2 ) To k B T 3 [exp(Haj/k B T Q ) - l] 3 

Thus, it is obtained that a > and a oc 1 /T approximately. Thus T c > T = (T L + T R )/2 is always 
true when the temperature difference between the left and the right thermal reservoirs is finite. 

Next, we will carry out numerical calculations for a three-terminal system shown in Fig. Q] 
The scalar model for the elastic wave is considered. And the model for thin geometry at low 



temperature is used so that the calculation is two-dimensional. So we can derive the transmission 



coefficient, r J( m , by the scattering matrix method 11261 . 
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In the calculation, we employ the 
following values of elastic stiffness constant and the mass density for GaAs J^oll : C44 = 5.99 x 10 1() 
Nm -2 and p = 53 17.6 kgm -3 , and choose W L = Wr = 10 nm, W T = 20 nm and W c = Dj = 10 nm. 
We truncate the sum of m in Eq. (OQ) at m = 10 Il29ll and limit the temperatures of the left and the 
right thermal reservoirs lower than T p h = finv/WtkB ~ 7.61 K (v = VcWp i s me sound velocity). 
At this low temperature, the phonon relaxation can be neglected I126n and the heat conduction is 
mainly determined by the ballistic transmission of the acoustic phonons. 

Fig. [2] shows T c - T vs AT for four different r s, where T is the averaging value of T L and T R , 
and T L = T + AT, T R = T - AT. First, it can be seen that when the difference between T L and T R 
is finite ( |2AT| > 0), the temperature of the central thermal reservoir T c is always higher than the 
averaging value of T L and T R , i.e. T c - T > 0, no matter which thermal reservoir has the higher 
temperature. Second, the temperature T c shows a quadratic dependence on AT, in agreement with 
Eq. ©. Third, as mentioned above that a oc I /To approximately, the curvatures of the curves 
depend strongly on the temperature T . The lower the temperature T , the larger the curvature. 

To study the nonlinear properties in the asymmetric three-terminal systems, we let the energy 
flux flows through the two terminals L and C and keep zero flux in terminal R (i.e. Q R = 0) and 
then calculate the temperature T R . Same as Eq. ©, with t cr m > T LR m , the temperature T R can be 
written as 



T R = T -/3(AT) + ^(AT) 2 + 0[(AT) 3 ], 



where 



7 = 



v f+°° ( d 2 n(T,a>) \ ,dco 



>o, 



> 0. 



(6) 



(7) 



(8) 



2m L R ,„ V^ffT 1 )^ n0J ( T CR,m + TLR,m)^ 

Fig. |3] shows T R - T vs AT for four different r s. The curves are also open up parabolic, in 
agreement with Eq. ©. The major difference with the symmetric case is that, here T R - T is no 
longer a symmetric function of AT with respect AT = but respect about AT = fi/y, and can 
be less than zero when < AT < 2/3/y. This means for some specific positive values of AT, the 
temperature T R can be less than the averaging value of T L and T c . The minimal value of T R is 
about T -/3 2 /2y. 



In summary, for symmetric three-terminal systems in the nonlinear response regime, we have 
found that the temperature of the central thermal reservoir is always higher than the averaging 
temperature of the left and the right thermal reservoirs. For the asymmetric three-terminal systems, 
the same nonlinear thermal properties can be also observed except that, in some special situations, 
the temperature of the central thermal reservoir can be lower than the averaging temperature of the 
left and the right thermal reservoirs. We would like emphasize that the temperature of the central 
thermal reservoir is insensitive to the details of the transmission characteristics. It is different 
from the electric case where the output voltage at the central branch shows fluctuations due to the 



transmission fluctuations[21]. This difference is due to the fact that the electric current is carried 



by a few electrons near the Fermi energy [31] but the heat flux is contributed by all phonons with 
different frequencies. The nonlinear thermal properties of the symmetric three-terminal systems 
can be used to control the heat flux. For example, when the temperature T c of the central thermal 
reservoir is lower than the averaging temperature (T L + T R )/2 of the left and the right thermal 
reservoirs, it is easier that the energy flux flows into the central thermal reservoir than the energy 
flux flows out from the central thermal reservoir when T c > (T L + T R )/2. 
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Figure Captions: 

Fig. Q] Schematic illustration of a symmetric three-terminal mesoscopic dielectric system. 

Fig. [2j T c - T vs A7\ calculated for a three-terminal system shown in Fig. [U for four different 
averaging temperatures, T = (T L + T R )/2, where T c , T L and T R are the temperatures of the 
central, the left and the right thermal reservoirs, respectively. All the temperatures are reduced by 
7> * 7.61 K. 

Fig. [31 T R - T vs AT, calculated for a three-terminal system shown in Fig.[H for four different 
temperatures, T = (T L + T c )/2, where T c , T L and T R are the temperatures of the central, the left 
and the right thermal reservoirs, respectively. All the temperatures are reduced by T ph « 7.61 K. 
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FIG. 2: 



8 




FIG. 3: 
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